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■^The  thrust  of  a  rocket  motor  Is  generally  calculated  from  the  stagnation 
pressure,  the  throat  area,  and  a  non-dimensional  thrust  coefficient  dependent 
upon  the  nozzle  expansion  are?  ratio.  A  comparison  is  made  between  momentum 
balance  and  pressure  integral  methods  of  calculating  the  ideal  vacuum  thrust 
coefficient,  and  the  results  used  to  examine  how  corrections  should  be  made  for 
nozzle  divergence,  skin  friction  and  two-phase  flow  to  obtain  the  real  thrust 
coefficient  of  the  motoi>^_l  _ .  -  ~~ 
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1  INTRODUCTION 

The  thrust  of  a  rocket  motor  is  generally  calculated  in  terms  of  the 
stagnation  pressure  (pQ)  and  the  throat  area  (Aa)  from  the  formula 

F  *  Cf  PQ  A* 

where  is  the  thrust  coefficient.  The  usual  procedure  is  to  calculate  the 

"ideal"  thrust  coefficient  from  one-dimensional  isentropic  flow  theory,  and  then 
to  apply  corrections  for  ambient  pressure  and  the  various  nozzle  "losses". 

The  "ideal"  thrust  coefficient  is  normally  calculated  by  considering  the 
overall  momentum  balance  of  the  motor.  It  may  also  be  calculated  by  considering 
the  pressure  distribution  over  the  internal  surfaces  of  the  motor.  In  this 
memorandum  the  two  methods  are  first  compared,  and  then  the  various  corrections 
are  identified  and  their  effect  on  the  final  value  of  the  thrust  coefficient  is 
cons i de red. 

2  IDEAL  THRUST  COEFFICIENT 

2.1  One-dimensional  isentropic  flow 

The  fundamental  equation  for  an  ideal,  one-dimensional,  compressible  flow 
is  Bernoulli's  equation  which  describes  the  conservation  of  energy: 


where  u  is  the  local  velocity,  a  is  the  speed  of  sound,  and  aQ  is  the 
speed  of  sound  at  stagnation  conditions.  If  we  remember  that 


a4  -  yRT 


then  this  can  be  transformed  to 


i  +^1h2 


(2.2) 


and  with  the  gas  law  and  isentropic  flow  relations 


-V 

p  ■  pRT  and  pp  ■  const 
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and 


p-2  ■  |>^"2] 


Y-1 


The  mass-flow  rate  at  any  section  is  constant: 


m  =  pAu  -  pA  Aa  ua 


where  a  *  subscript  denotes  the  throat  conditions 


Ma  -  1 


(2.3) 


(2.4) 


(2.5) 


The  mass-flow  rate  can  therefore  be  calculated  by  substituting  (2.1)  and  (2.3) 
in  (2.5)  to  give 


It! 


(2.6) 


Finally,  we  can  use  (2.5)  to  give  relationship  between  Mach  number  and  area  ratio 


A_ 

A* 


(2.7) 


2.2  Momentum  balance  equation 

If  we  consider  the  overall  momentum  balance  of  the  control  volume  shown 
in  Fig.  1: 


-  A  u  +  p  A 


F 
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where  u  is  the  gas  velocity  at  exit,  p  the  static  gas  pressure  and  A 
6  e  e 

the  exit  cross-sectional  area.  For  the  moment  we  shall  assume  that  the  ambient 
pressure  is  zero. 

From  (2.1) 


u 

e 

a 

o 


2 


o 


(2.8) 


where  Me  denotes  the  Mach  number  at  the  exit  plane. 


Also,  combining  (2.1*)  and  (2.7): 


2  ^2TFT7 

.  yTT  j 


-i 


(2.9) 


so  that,  using  (2.6),  (2.8)  and  (2.9) 


F 


} 


or 


(2.10) 


2.3  Pressure  integral  method 

While  the  method  shown  above  provides  a  simple  means  of  calculating  the 
ideal  thrust  coefficient,  it  leaves  some  ambiguities,  particularly  when  we  come 
to  consider  how  various  "losses"  In  the  thrust  coefficient  arise.  Most 
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immediately  we  observe  that  the  final  term  (pg  Ag)  in  the  momentum  equation 
does  not  represent  a  force  acting  on  the  physical  structure  of  the  nozzle  at 
all,  but  is  really  an  "adjustment"  made  to  the  equation  because  the  control 
volume  has  been  drawn  across  the  exit  plane  of  the  nozzle.  The  influence  of 
the  nozzle  will  actually  extend  downstream  of  this  point  within  an  envelope 
bounded  by  the  Mach  lines  extending  from  the  lip  of  the  nozzle  to  the  centre 
line.  In  part  this  ambiguity  arises  because  the  "one-dimensional  flow"  assump¬ 
tion  cannot  be  true  if  the  flow  is  to  expand. 

The  alternative  approach  is  to  observe  that  the  thrust  is  generated  by  the 
integral  of  the  internal  pressure  over  the  surface  of  the  nozzle: 

F  =  |  pdA 

where  A  is  the  axial  cross-section  of  the  flow.  In  fact  the  integration  will 
have  to  be  made  in  two  parts  -  one  over  the  head-end  of  the  motor,  and  the 
second  over  the  length  of  the  nozzle.  At  this  point  we  can  also  allow  for  the 
finite  contraction  ratio  of  the  nozzle 


F  =  p  A  +  pdA 

00  Ja 


(2.11) 


If  we  differentiate  the  mass-flow  equation  (2.9)  we  have 


+  dA  +  du  _  Q 
p  A  u 


(2.12) 


and  since  the  differential  forms  of  the  gas  law  and  isentropic  flow  equations 
are 

d£  .  d£  .  dT  _  0 
p  P  T 

and 


d £ 

P 


d£ 


dp  1  dT 

P  "  TFTT  “ 
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which,  by  differentiating  (2.1)  becomes 


dp 

P 


u  da 


a 


M 


2  d_u 
u 


Equation  (2.12)  can  therefore  be  written  as 


£  -  <H2  -  1)  ^ 

A  u 


(2.13) 


which  also  indicates  how  the  nozzle  must  change  from  a  contraction  to  an 
expansion  as  the  Mach  number  passes  unity  to  maintain  a  positive  acceleration. 

In  fact  we  can  use  equation  (2.1)  to  convert  (2.13)  entirely  to  terms  of 
Mach  number: 

dM  _  du  _  da 
M  u  a 

and 

2 

u  du  +  — r  •  a  da  ■*  0 

y-1 


and 


dA  M2  -  1  dM 

‘  " 


(2.1A) 


which  is  the  differential  form  of  (2.7)*  We  can  now  substitute  this  equation 
together  with  (2.7)  and  (2.1*)  in  (2.11): 


F 


P  A 
Ko  o 


Po  A* 


2 

Y+1 


1+1 

2 (y-1) 


rM 
I  e 


M 

o 


Lstzll 


dM 
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o 

The  ratio  Aq/A.v  will  be  determined  by  the  entry  Mach  number  to  the  nozzle, 
as  in  equation  (2.7).  When  the  entry  Mach  number  is  small,  this  term  cancels 
with  the  first  term  in  the  integral,  as  might  be  expected  since  C^.  is  not 
very  sensitive  to  entry  Mach  number.  The  integral  then  proves  to  be  identical 
with  (2.10) 


A  f 

m  *T+  ( 


2 

Y+1 


1 2 (y- 1 ) 


C 


f 


Y+1 

2Ty-TT 


(1  +  Y  Me2) 


3  CORRECTION  FOR  AMBIENT  PRESSURE 

So  far  we  have  assumed  that  the  ambient  pressure  is  zero.  In  practice 
with  a  finite  ambient  pressure  (p  )  we  should  continue  the  integral  in 

a 

section  2.3  over  the  outside  of  the  motor.  The  integral  is  simplified  by  the 
fact  that  p  may  be  assumed  everywhere  constant.  The  thrust  is  then  given  by 

Cl 


T  =  Cf  PQ  A*  + 


r° 


pa  dA 


Ce  p„A*  -  p  A„ 
f  ro  *  ra  e 


The  integral  here  goes  from  the  lip  of  the  exit  cone  to  the  centre  of  the  head 
end  of  the  motor.  Thus 


-  C 


fv 


Pa  Ae 
Po  A* 
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If  p  is  not  constant  over  the  outside  of  the  motor,  e.g.  due  to  induced  flow 
past  the  motor,  then  it  will  be  necessary  to  re-evaluate  the  integral  to  allow 
for  this  effect- 

14  THRUST  COEFFICIENT  LOSSES 

4. 1  Genera  1 

There  are  a  number  of  losses  to  the  overall  thrust  which  do  not  affect  the 
thrust  coefficient,  and  we  shall  not  consider  those  here.  Under  this  category 
we  shall  ignore  effects  which  alter  the  flow  area  of  the  throat  or  the  upstream 
stagnation  conditions.  The  major  losses  to  thrust  coefficient  are  then 

-  divergence  losses 

-  skin  friction  losses 

-  two-phase  flow  losses 

We  must  now  consider  how  these  losses  are  to  be  integrated  into  calculation  of 
the  thrust  coefficient. 

4.2  Divergence  losses 

On  exit  from  the  nozzle  the  flow  has  a  component  of  momentum  in  the  radial 
as  well  as  the  axial  direction.  Malina  £lj  calculated  the  loss  of  thrust  to 
radial  momentum  by  assuming  that  the  exit  flow  formed  a  cone  of  streamlines 
originated  from  a  point  source  near  the  throat.  He  showed  that  the  reduction 
in  axial  momentum  in  this  instance  was  by  a  factor 

1  +  cose 
nd  2 

where  0  is  the  half  angle  of  the  cone.  Nowadays  more  complete  computer  flow 
predictions  are  available,  but  it  appears  that  Malina's  formula  provides  a  high 
degree  of  accuracy  if  the  half-angle  of  the  flow  at  exit  is  used. 

Using  Malina's  formula  one  might  assume  that  only  the  velocity  terms  in  the 
momentum  balance  equation  would  be  affected.  However,  Landsbaum  ["  2 J  pointed 
out  that  the  exit  area  involved  in  Malina's  assumption  is  actually  the  cap  of  a 
sphere  and  not  the  flat  exit  plane,  and  that  the  "one  dimensional"  calculation 
of  Cf  is  initially  high,  and  using  Malina's  correction  on  all  the  terms  in 
C^v  provides  a  closer  approximation  to  the  actual  correction,  the  errror  then 
being  about  0.1$  pessimistic  for  a  10:1  area  ratio  nozzle  with  0  »  15°. 
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k.},  Skin  friction  loss 


The  drag  force  acting  on  a  nozzle  expansion  cone  can  be  calculated  as: 

f  C  r  O 


Drag 


p  u  2n  rdx 


where  r  is  the  local 
coefficient.  The  skin 
boundary  layer  methods 
as  a  first  approximati 


radius  of  the  nozzle,  and  c^.  is  a  skin  friction 
friction  coefficient  can  be  calculated  by  compress 
which  include  the  effects  of  the  accelerating  flow 
on 


i  ble 
,  but 


^  =  0.0256  Re  ' ‘2 

2  x 


where  Re^  is  the  Reynolds  number  based  on  distance  down  the  flow.  We  can 
observe  that  c^  decreases  as  we  move  further  downstream. 

I f  we  observe  that 


rti  =  pu  •  ir  r 


Pn  A* 


Drag 


po  A* 


C^ 


lidx 

r 


Now  u/c.,.  is  a  component  of  the  thrust  coefficient  momentum  balance,  and  like 
the  thrust  coefficient  it  increases  only  slowly  as  we  go  down  the  nozzle,  with 
a  value  of  about  0.7  at  the  throat  increasing  to  about  1.65  for  a  12:1  area 
ratio  nozzle.  The  increase  in  u/c*  is  actually  rather  greater  than  the 
decrease  in  skin  friction  coefficient,  but  it  is  probably  close  enough  to  write 


Drag 
Po  A* 


=  c 


—  f  - 

fe  c*  j  r 


dx 


where  f  is  a  non-dimensional  factor  accounting  for  the  variation  of  c^  u/c* 
with  x  .  We  therefore  see  that  if  we  define  a  skin  friction  loss  efficiency 


1  -  c 


fe 


f  f!* 
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it  will  apply  only  to  the  velocity  term  in  the  thrust  coefficient.  Typically 
rif  is  about  0.99* 

4 . 4  Two-phase  flow  losses 

Most  particularly  in  the  case  of  aluminized  composite  propellants  (where 
up  to  25%  of  the  flow  can  be  in  the  form  of  A^O^  particles)  the  presence  of 
solid  particles  in  the  flow  induce  flow  losses  both  by  retaining  energy  which 
should  contribute  to  the  thermodynamic  expansion  ("thermal  lag")  and  by  having 
a  lower  acceleration  than  the  gas  ("velocity  lag").  If  a  fraction  £  of  the 
exhaust  is  solid  particles,  the  thrust  is 

Thrust  =  (1  -  £)  m  ug  +  pe  Ae  +  S  ™  up 

where  u^  is  the  velocity  of  particles  at  the  exit  plane.  If  the  velocity 
lag  is  £ 


T  =  (1  -  £  +  £c)  m  ue  +  pe  Ag 

This  simplified  analysis  suggests  that  the  correction  for  two-phase  flow 
should  only  be  applied  to  the  velocity  term.  However,  it  is  a  simplified 
analysis.  If  we  examine  in  more  detail  what  is  happening  within  the  nozzle 
we  find  that  the  drag  of  the  solid  particles  will  slow  the  velocity  of  the  gas 
phase  and  increase  its  temperature,  thereby  reducing  the  Mach  number  of  the  gas 
The  thermal  lag  of  the  particles  will  also  tend  to  add  heat  to  the  flow  down¬ 
stream.  Since  in  the  expansion  cone  at  least  thrust  is  generated  only  by  the 
pressure  of  the  gas,  a  reduction  by  Mach  number  represents  a  reduction  in  the 

whole  of  C, 
fv 

Part  of  the  problem  in  assigning  a  correction  for  two-phase  flow  arises 
from  the  problems  of  calculating  the  effect  accurately.  A  large  part  of  the 
loss  will  occur  upstream  of  the  nozzle  throat,  and  so  will  affect  the  character 
istic  velocity  c  also.  It  is  significant  that  where  exact  formulations  of 
two-phase  flow  problems  have  been  carried  out  (e.g.  Marble  £  3j)  these  have 
concentrated  on  calculating  the  reduction  in  specific  impulse,  combining  C^ 
and 

5  CONCLUSIONS 

This  memo  has  shown  that  the  value  for  the  ideal,  "one-dimensional"  vacuum 
thrust  coefficient  is  the  same  when  calculated  by  a  momentum  balance  or  by  a 
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pressure  distribution  method.  Investigation  of  the  various  losses  suggests 
that  the  real  thrust  coefficient  can  be  calculated  from: 


nd  nt  [nf 


P  A 

c,  +  0  -  n,)  e-.-e 
fv  f  Pn  A* 


] 


A 

e 

\ 


whe  re 


is  the  divergence  loss  efficiency 
nt  is  the  two-phase  flow  loss  efficiency,  and 
n^.  is  the  skin-friction  loss  efficiency. 
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